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Abstract 

Pulsar timing experiments aimed at the detection of gravitational radiation have been performed 
for decades now. With the forthcoming construction of large arrays capable of tracking multiple 
millisecond pulsars, it is very likely we will be able to make the first detection of gravitational 
radiation in the nano-Hertz band, and test Einstein's theory of relativity by measuring the 
polarization components of the detected signals. Since a gravitational wave predicted by the 
most general relativistic metric theory of gravity accounts for six polarization modes (the usual 
two Einstein's tensor polarizations as well as two vector and two scalar wave components), we 
have estimated the single- antenna sensitivities to these six polarizations. We find pulsar timing 
experiments to be significantly more sensitive, over their entire observational frequency band 
(~ 10^^ — 10^^ Hz), to scalar-longitudinal and vector waves than to scalar-transverse and tensor 
waves. At 10~^ Hz and with pulsars at a distance of 1 kpc, for instance, we estimate an average 
sensitivity to scalar-longitudinal waves that is more than two orders of magnitude better than the 
sensitivity to tensor waves. Our results imply that a direct detection of gravitational radiation by 
pulsar timing will result into a test of the theory of general relativity that is more stringent than 
that based on monitoring the decay of the orbital period of a binary system. 

PACS numbers: 98.80.-k, 95.36. +x, 95.30.Sf 
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I. INTRODUCTION 



The pulsar timing technique for gravitational wave searches relies on the pulsar emission 
of highly periodic radio pulses, which are received at the Earth and cross-correlated against 
a template of the pulsed waveform. A gravitational wave (GW) propagating across the 
pulsar-Earth radio link introduces a time shift in the received pulses that is proportional to 
the amplitude of the GW. Since pulsars, and millisecond pulsars in particular, are the most 
stable clocks in the Universe over time-scales of years, they provide a unique element for 
performing gravitational wave searches in the nano-Hertz band. 

Pulsar timing sensitivities have improved significantly over the years, and with the advent 
of the forthcoming arraying projects it is very likely we will be able to make the first detection 
of gravitational radiation in the nano-Hertz band. The first unambiguous detection of a 
gravitational wave signal will also allow us to test Einstein's general theory of relativity by 
measuring the polarization components of the detected signals [HE]- Since general relativity 
is the most restrictive among all the proposed relativistic metric theories of gravity [3J, as 
it allows for only two of possible six different polarizations P], by asserting that the spin-2 
("tensor") polarizations are the only polarization components observed, we would make a 
powerful proof of the validity of Einstein's theory of relativity. Corroboration of polarization 
measurements with estimates of the propagation speed of the observed gravitational wave 
signal will provide further insight into the nature of the observed radiation and result into 
the determination of the mass of the graviton. 

The advent of the "dark energy and dark matter" problem, the proposal of new 
inflationary scenarios, the theoretical attempts to quantize gravity, and the possible existence 
of extra dimensions [3], have all stimulated an increasing interest in theories of gravity that 
are alternatives to Einstein's theory of general relativity. Among the proposed theories, 
scalar-tensor theories predict GWs with one or two polarization modes with helicity s = 
besides the two usual modes with helicity s = ±2 [H [3]. Also theories that introduce 
geometrical corrections in the Einstein-Hilbert Lagrangian have appeared in the literature. 
The so-called f{R) theories are of this kind and predict the existence of two additional scalar 
degrees of freedom pj. Vector polarization modes, on the other hand, can appear in the 
"quadratic gravity" formulations [6], and in the context of theories for which the graviton 
has a finite mass (such as the Visser theory [7j). In fact, this last group of theories are the 
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most general in that they allow for GWs with six polarization modes. 

The problem of observing additional polarization modes with arrays of pulsar timing 
has been analyzed by Lee et al. [2j in the context of searches for an isotropic stochastic 
background of gravitational radiation. Their approach relied on cross-correlating pairs 
of time of arrival residuals (TOAR) from an array of pulsars. In what follows we will 
focus instead on the single-antenna pulsar timing sensitivity to GWs characterized by the 
presence of all possible polarization modes. This work is a follow up to our recent derivation 
of the Laser Interferometer Space Antenna (LISA) sensitivities to GWs predicted by the 
most general relativistic metric theory of gravity [8]. There we showed, in particular, 
that the LISA sensitivity to scalar-longitudinal waves is significantly better than that to 
tensor waves when the wavelength of the gravitational wave signal is shorter then the size 
of the detector. This is because the root-mean-squared of the one-way Doppler responses to 
waves with longitudinal components experience an amplification that is roughly proportional 
to the product between the frequency of the observed signal and the distance between 
the spacecraft. In the case of pulsar timing experiments the resulting amplification to 
waves with longitudinal components over those purely transverse is rather significant due to 
their operational frequency bandwidth (10^^ — 10^^ Hz) and typical pulsar-Earth distance 
1 kpc) [2j. In the same spirit of 0, here we evaluate the pulsar timing sensitivities 
following the standard definition used for all other GW detectors: the strength of a sinusoidal 
gravitational wave signal, averaged over the sky and polarization states, required to achieve 
a given signal-to-noise ratio (SNR) over a specified integration time, as a function of Fourier 
frequency. 

The paper is organized as follows. In Sec. |TT] we provide the expression for the TOAR 
response to GWs with helicity s = 0, ±1 and ±2. Since the time-derivative of the TOAR 
is equal to the Doppler response, we do not derive its expression here as this was presented 



in our earlier publication [8]. In Sec. Ill, after summarizing the main noise sources affecting 
pulsar timing experiments, we evaluate the sensitivities of these experiments to each of the 
different GW polarizations. We find pulsar timing experiments to be two to three orders 
of magnitude more sensitive, over their entire observational frequency band (~ 10"^ — 10~^ 
Hz), to scalar-longitudinal waves than to tensor waves. This result implies that a direct 
detection of gravitational radiation by pulsar timing will allow us to perform a dynamical 
test of the theory of general relativity that is more stringent than that based on monitoring 



the decay of the orbital period of a binary system (see for instance [9]). Finally in Sec. IV 
we summarize our results and present out conclusions. 

Throughout the paper we will be using units such that c = 1 except where mentioned 
otherwise. 

II. THE PULSAR TIMING RESPONSE 

Since the time-derivative of the TOAR is equal to the relative frequency change of the 
radio pulses, in what follows we derive the pulsar timing sensitivity by relying on the one- 
way Doppler responses to the noises and to a GW signal containing all the six polarizations. 
Although the one-way Doppler response to a general GW was first derived by Hellings [10] , a 
simpler and more compact expression for it was recently obtained by us which is identical 
in form to that first obtained by Wahlquist [TT] for GWs with s = ±2. 

Our notation is summarized by Figjl} hL is the vector oriented from Earth to the pulsar, 
which continuously emits radio pulses that are received at the ground station; k is the unit 
wave- vector of a GW propagating in the +z direction, and 6 and are the usual polar angles 
associated with k. 

If we denote with ycwit) the relative frequency changes induced by a gravitational wave 
signal on the radio link between a pulsar and the Earth, its expression is equal to [S] [12] 

VGwit) = {l-k-n) [^{t -{1 + k- n)L) - ^(t)] , (1) 

where 

^(t) = — =r . (2) 

2 l-{k-n)^ 

In the above expression hij{t) are the spatial components of the GW metric perturbation 
corresponding to the following space-time line element 

ds'^ = -dt^ + {5ij + hijit - z))dx'dx^ , (3) 

and \hij\ ^ 1. In this way, the only restriction on h^y is that its temporal components 
are null, i.e., h^^ = and therefore it contains six degrees of freedom. A general GW 
perturbation can be written as a sum of six components in the following way 
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h,^{t-z) = J24^h^r){t-z) , (4) 

r=l 



z 




FIG. 1. The radio pulses emitted by the pulsar are received at Earth by a radio telescope. The 
gravitational wave train propagates along the z direction, and the two polar angles {0, (p) describe 
the direction of propagation of the radio pulses relative to the wave. See text for a complete 
description. 



where e\j are the six polarization tensors associated with the six waveforms of the 
gravitational wave signal. In the Cartesian coordinates system described in Fig. [T| the 
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above polarization tensors assume the following matricial forms 

^ l\ 




'^O 0^ 




1 



-(2) 



1 



'(3) 



'^O 0^ 

1 
1 



= (4) 
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'^O 1 0^ 
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'^l 0^ 
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By analyzing the behavior of the above matrices under rotation around the z axis we find 



that the polarization tensors e[]^ and e^"'' have helicity s = and we call them scalar- 



.(6) 



longitudinal and scalar-transversal modes respectively; e[^^ and e\j' have helicity s 



,(3) 



±1 



and they are called vector-modes. Finally the usual tensor modes are represented by the 
polarization tensors e\f and e^f, which have helicity s = ±2. 



III. SENSITIVITY 

The sensitivity of a gravitational wave detector has been traditionally taken to be equal to 
(on average over the sky and polarization states) the strength of a sinusoidal gravitational 
wave required to achieve a given signal-to-noise ratio (SNR) over a specified integration 
time, as a function of Fourier frequency. Here we will assume a SNR = 1 over an 
integration time of 10 years. Explicitly, we have computed the sensitivity using the following 
formula: ^y[Sy{f)B]/ (rms of the GW response), where Sy{f) is the spectrum of the relative 
frequency fluctuations due to the noises affecting the pulsar timing data, / is the Fourier 
frequency, and B is the bandwidth corresponding to an integration time of 10 years. The 
expression for the spectrum of the noise Sy{f) we used in our sensitivity calculations relies on 
the noise-model discussed in [13]. In particular, we have assumed that: (i) multiple-frequency 
measurements can be implemented in order to adequately calibrate timing fluctuations due 
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to the intergalactic and interplanetary plasma, and (ii) the tracked millisecond pulsars have 
frequency stabilities better than those of the operational ground clocks. A recent stability 
analysis of presently known millisecond pulsars [2] has shown that there might exist some 
with frequency stabilities superior to those displayed by the most stable operational clocks 
in the (10~^ — 10~^Hz) frequency band. 

Under these assumptions, the expression for the noise spectrum, Sy{f), is characterized 
by the contribution due the ground clock and a white-timing noise of 100 nsec. in a Fourier 
band +/- 0.5 cycles/day (i.e. one sample per day). The 100 nsec. level is the current timing 
goal of leading timing array experiments as three pulsars are being timed to this level [H] . 
Following Jenet et al. pjj the expression for the noise spectrum Sy{f) is equal to 

Sy{f) = [4.0 X 10-^^ f-^ + 3.41 X 10-^ f]Hz-^ . (6) 

We have assumed the vector waves to be elliptically polarized and monochromatic, with their 
waveforms, {h^'^\ h^^^), written in terms of a nominal amplitude, H, and the two Poincare 
parameters, ($,r), in the following way 

h^^\t) = Hsm{T) sm{ut + $) , 

h'-^^t) = H cos{T) sm{ujt) . (7) 

The sensitivity for elliptically polarized tensor components can be written in the same way 
as the vector modes by just replacing {h^'^\ h^^^) with {h^'^\ h^^'^) in Eqs. For scalar 

signals instead, the two wave functions, /i*^^)), have been treated as independent since 

one is purely longitudinal and the other purely transverse to the direction of propagation of 
the GW. 

Since we are interested in the sensitivity averaged over the sky, we calculated the averaged 
modulus squared of the Fourier transform of the GW response taken over the sky by assuming 
sources uniformly distributed over the celestial sphere; in the case of tensor and vector signals 
we also averaged over elliptical polarization states uniformly distributed on the Poincare 
sphere for each source direction. Fig. |2] shows the GW sensitivities corresponding to the 
noise spectrum given in Eq. ^ and after assuming a coherent integration of ten years 
{B = 1 cycle/10 years). The sensitivities of pulsar timing to vector and scalar- longitudinal 
polarized waves are significantly better than to tensor and scalar-transverse waves over the 
entire observational frequency band (10^^ — 10^^ Hz). As explained in our earlier publication 
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Pulsar Timing Sensitivities 




iog(f(Hz)) 



FIG. 2. Sensitivities of pulsar timing experiments to gravitational waves with tensor {s = ±2), 
vector (s = ±1), and scalar (s = 0) components. We may notice how more sensitive pulsar timing 
experiments are to vector and scalar-longitudinal signals than to tensor and scalar-transverse waves. 
This is true over the overall accessible frequency band, and more pronounced at higher frequencies. 
See text for more details. 

[8] , where this effect was first noticed in the context of the LISA mission, the physical reason 
behind this sensitivity enhancement is due to the difference in the amount of time a "pulse" 
tensor wave and pulse-wave with a longitudinal component affect the pulsar-Earth radio 
link. A tensor signal propagating orthogonally to the pulsar-Earth radio link (direction for 
which the one-way Doppler response can reach its maximum magnitude in this case) will 
only interact with the electromagnetic link for the time it takes the wave front to cross the 
radio link. On the other hand, if (for instance) a scalar-longitudinal wave propagates along 
the pulsar-Earth direction (over which the Doppler response will achieve its maximum in 
this case), the frequency of the radio link will be affected by the gravitational wave for the 
entire time L it takes the wave to propagate from the pulsar to Earth, resulting into an 
amplification of the Doppler response. 
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In mathematical terms the argument can be explained in the following way. For the 
tensor modes, the maximum response occurs when k ■ h ^ and it is proportional to 
|sin(7r/L)|. For the scalar-longitudinal mode instead the maximum response is achieved 
when A; ■ n — )■ 1 and is proportional to fL [15] . When the average over all possible directions 
and polarizations of the incoming waves is performed, the pulsar timing sensitivity to scalar- 
longitudinal waves is better than that to tensor waves by more than one order of magnitude 
at / = 10~^Hz and three orders of magnitude at / = 10~^Hz if we assume a pulsar out to a 
distance of 1 kpc. The same explanation applies to vector waves since they also affect space 
along their direction of propagation. 

Since the measurements of the decay of the orbital period of the binary pulsars PSR 
B1913+16 and PSR B1534+12 were performed at an accuracy level of about 1 percent by 
comparing them against the expected energy loss due gravitational wave emission predicted 
by Einstein's theory of relativity, we infer that these measurements resulted into a 10 percent 
upper limit assessment of the amplitudes of the non-Einsteinian polarizations. This is 
because the energy of a GW signal is proportional to the square of its amplitude and we have 
also assumed that the mentioned accuracy can be related to the ability of this observations 
to access the non-tensor polarizations. In this sense, the direct detection of a gravitational 
wave signal with pulsar timing could potentially achieve a much higher level of accuracy in 
testing the theory of relativity because of its enhanced sensitivity to non-Einsteinian modes. 
As an example, consider an hypothetic detection of a gravitational wave signal around the 
Fourier frequency 10~^ Hz. If an appropriate data analysis indicates that only the tensor 
components are present in the data, from our estimated sensitivities we will conclude that, 
on the average, the amplitude of a scalar longitudinal wave must be smaller than that of 
the tensor wave by a factor of about 300 or more. This translates into a test of the theory 
of relativity that is better than 0.3 percent in wave's amplitudes. Detections occurring at 
higher frequencies would result into more stringent tests of the theory, as implied by Fig. |2} 

As a final comment, our estimated sensitivity enhancements experienced by the vector and 
scalar-longitudinal waves over tensor waves are independent of the assumptions underlining 
the spectrum of the noise used in our calculations. 
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IV. CONCLUSION 



The main result of our work has been that pulsar timing experiments are more sensitive 
to scalar-longitudinal and vector signals than to scalar-transverse and tensor waves over 
their overall accessible frequency band. In particular, the fact that the sensitivities of 
pulsar timing experiments to waves with longitudinal components are better than to purely 
transverse signals by two to three orders of magnitude indicates that future experiments 
performed by forthcoming arraying projects will be able to assess the polarization of the 
detected waves. This will result into a dynamical test of Einstein's theory of relativity 
significantly more stringent than those performed to date based on measuring the decay of 
the orbital period of a binary system containing a pulsar. 

The ability of entirely reconstructing the waveforms of a detected signal will require the 
use of six or more millisecond pulsars [10], as it can be argued from a simple counting 
argument of the number of unknowns characterizing the most general gravitational wave 
signal. This task will be most suitably addressed by forthcoming pulsar arraying projects, 
and it will be the topic of our future investigation. 
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